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Abstract 

The two-dimensional Schrodinger operator with a uniform magnetic 
field and a periodic zero-range potential is considered. For weak mag- 
netic fields and a weak coupling we reduce the spectral problem to the 
semiclassical analysis of one-dimensional Harper-like operators. This 
shows the existence of parts of Cantor structure in the spectrum for 
special values of the magnetic flux. 
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1 Introduction 



The spectral properties of a charged particle in a two-dimensional system 
submitted to a periodic electric potential and a uniform magnetic field cru- 
cially depend on the arithmetic properties of the number 9 representing the 
magnetic flux quanta through the elementary cell of periods, see e.g. [8] for 
a description of various models. Since the works by Azbel [7J and Hofstadter 
[2U] it is generally believed that for irrational 9 the spectrum is a Cantor 
set, and the graphical presentation of the dependence of the spectrum on 9 
shows a fractal behavior known as the Hofstadter butterfly. (We note that 
fractal spectral diagrams has been found recently in the three-dimensional 
situation as well, see e.g. |9J). After intensive efforts during more than 
twenty years this was rigorously proved recently (Ten Martini conjecture) for 
all irrational values of 9 for the discrete Hofstadter model, i.e. the discrete 
magnetic Laplacian admitting a reduction to the almost Mathieu equation, 
see [3] and references therein. 

Only few results are available for other models. Traditionally, semiclassical 
methods have played an important role in the analysis of the two-dimensional 
magnetic Schrodinger operators with periodic potentials, see e.g. [TO] for a 
review. In particular, the bottom part of the spectrum for strong magnetic 
fields can be described up to some extent using the tunneling asymptotics 
in a very general setting [21]. Concerning a more detailed analysis, it was 
shown by Helffer and Sjostrand [251 E6J, [27] that the study of some parts 
of the spectrum for the Schrodinger operator with a magnetic field and a 
periodic electric potentials reduces to the spectral problem for an operator 
pencil of one-dimensional quasiperiodic pseudodifferential operators (see be- 
low, Section [3]); recently this correspondence was extended up to a unitary 
equivalence [I5[ [16]. Under some symmetry conditions for the electric po- 
tentials, the operator pencil reduces to the study of small perturbation of 
the continuous analog of the almost-Mathieu operator, which allowed one 
to carry out a rather detailed iterative analysis for special values of 9. In 
particular, in several asymptotic regimes a Cantor structure of some part of 
the spectrum was proved. 

In the present paper we are interested in the spectrum of the two-dimensional 
magnetic Schrodinger operator with periodic zero-range potentials (called 
also point perturbations). Perturbations of such a kind became a rather 
popular model in quantum mechanics due to a possibility of an analytic 
investigation whose results are in a good agreement with experiments, see 
the monographs [21 [13]. Various aspects of the spectral analysis of magnetic 
Hamiltonians with point perturbations were discussed in numerous works, 
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see e.g. [H El HU El HE1 ESI [221 E2] and references therein. Our principal 
aim in this paper is to show that the technique of Helffer and Sjostrand 
is still applicable and is general enough to handle this non-standard class 
of operators; such a correspondence between the semiclassical analysis and 
solvable models seems to appear for the first time. This provides a proof 
of the existence of Cantor parts in the spectrum in a certain asymptotic 
situation. 

Our starting point will be the construction of [IB] for the resolvent of the per- 
turbed operator, which reduces the spectral problem to the study of a family 
of discrete operators (see Section [2]). We then construct a family of pseu- 
do differential operators (effective Hamiltonians) having the same spectrum 
as the discrete operators (Section HJ). In Section [5] we obtain some estimates 
for the Green function of the Landau Hamiltonians, which are used to show 
that in the weak magnetic field limit the effective Hamiltonians obey the 
conditions needed for the Helffer-Sjostrand analysis (Section [6] and the main 
result in Theorem 16 Ah . 

2 Magnetic operator with zero-range poten- 
tials 

We start with the two-dimensional Schrodinger operator with a uniform mag- 
netic field in the Landau gauge (Landau Hamiltonian) , 



Recall that the spectrum of consists of the infinitely degenerate eigenval- 
ues E n = (2n — l)h, n e N, called the Landau levels. We are going to study 
periodic zero-range perturbations of Hh supported by the set 1? . Physically, 
such operators model periodic arrays of identical small impurities. The in- 
teractions can be intuitively understood as operators corresponding to the 
limit for e — > of the operators 



for localized potentials V e approaching the Dirac 5 function and a special 
choice of the constants a(e). An alternative (but essentially equivalent) ap- 
proach, which is more suitable for the spectral analysis, involves the use of 
self-adjoint extensions [12J. In this case one considers first the restriction of 




h>Q. 



H(s)=H h + a(e) £ V e (--m), 



meZ 2 
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Hh to the functions vanishing at all points of Z 2 (this operator is well defined 
as all functions in the domain of H h are continuous), to be denoted by S. By 
zero-range perturbations of Hh, one means then self-adjoint extensions of S. 

We restrict ourselves by considering the standard one-parameter family of 
operators (self-adjoint extensions) Hh, a , «eK \ {0}, whose domain consists 
of the functions / having logarithmic singularities at the points of Z 2 , 

f(x) = ^ log + b m (f) + o(l) as \x-m\= o(l), 

Z7T \x — ml 
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a m (f),b m (f) G C, m G Z : 

and satisfying the "boundary conditions" a m (f) + ab m (f) = 0, m G Z 2 . 
The action of the operators in the sense of distributions is given by the 
following expression (Fermi pseudopotential): 

H h ,af(x) := H h f(x) + a ^ S(x-m)(l + log \x - m\ (x - m, V))/(x), 

where 5 is the Dirac delta function. The above boundary conditions sat- 
isfy the natural assumption of form-locality [21], i.e. the associated bilinear 
form vanishes on functions with disjoint supports. This helps to avoid some 
pathological spectral properties appearing for more general boundary condi- 
tions EDI. 



Let us repeat the main constructions of [H] to show how to handle the 
spectral problem for Hh, a - Note first that the Green function (the integral 
kernel of the resolvent) of is known explicitly, 

1 / ih(x 1 -yx)(x 2 +y2) h(x-yf 



G h (x,y; z) = — exp ( - 



4tt 1 V 2 



where U is the Kummer confluent hypergeometric function. The Green func- 
tion plays a crucial role in the spectral analysis of the operators Hh, a - Namely, 
let us define, for z spec Hh, 

1 . 1 



q(z, h) := lim (G h (x, y\ z) - — log 



\x-y\ 



1 / , /l z n , h , A 



4vrV rv 2 2h J °2 

where ip is the logarithmic derivative of the T-function. Note that the limit 
is independent of the choice of y G M 2 . For z ^ spec Hh define an operator 
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Q(z, h) : £ 2 (Z 2 ) — > £ 2 (Z 2 ) given in the canonical basis by the matrix 

jq(z,h), m = n, 

Q{m, n; z,n) = < 

I Gh(m, n; z), otherwise. 

Then the Green function G^, a of ff^Q, is given, for z £ specify U specify, 
by 

G hta (x,y;z) = G h (x,y;z) 

-a 7] (a<2(z, /i) + 1) 1 {m,n)G h (x,m]z)G h (n,y;z), (2.3) 

where the series on the right-hand side converges in the strong resolvent 
sense. An important ingredient of the above formula (12.31) is the relation 

spec H h a \ spec H h = {2 ^ spec H h : G spec (Q(z, h) + aT 1 ) } . (2.4) 

As the set spec H h is discrete (the set of the Landau levels), Eq. (12 .4p provides 
an almost complete characterization of the spectrum of Hh, a , and this is the 
starting point for the subsequent analysis. Of crucial importance will be also 
the identity 

Q(m + k,n + k;z,h)=e- ihk2{mi - ni) Q(m,n;z,h), m,n,keZ 2 , (2.5) 

which can be verified directly; this expresses the invariance of Hh under the 
magnetic translations [34J. 

3 Strong type I operators 

In the works [25j [26j [27] a machinery was developed for an iterative semi- 
classical analysis of a special class of pseudodifferential operators. One was 
concerned with the non-linear spectral problem (or, in other words, with 
the spectral problem for an operator pencil). Namely, for a family of self- 
adjoint operators A(fi) depending on a real parameter /i, one means by the 
/x-spectrum //-spec A(p) the set of all \x such that G spec A(p). In particu- 
lar, the //-spectrum of the family A(f/,) := A — /i for a self-adjoint operator 
A is exactly the spectrum of A. 

Note that in the above terms the problem (12.41) is exactly to find the z- 
spectrum for the family Q(z, h) + a -1 . 
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Consider a bounded function (symbol) L : R 2 — > C, L = L(x,p). By the 
Weyl quantization procedure one can assign to L an operator L h in £ 2 (M), 

lkf{x) = ^ LI/ p{x ~ y)/hL (^' p)f{y)dpdy - (3 - i} 

Note that L can depend itself on h and other parameters as well. The 
operator Lh obtained is referred to as the Weyl /i-quantization of L (or the 
/i-pseudodifferential operator with the symbol L), and quantum Hamiltonians 
resulting from symbols periodic in both x and p are often called Harper-like 
operators. In particular, the symbol L(x,p) := cosp + cosx produces the 
Harper operator on the real line, 

£*/M = /(* + *>+/('-*> + cos , /(*). (3.2) 

In [27], in order to treat the Harper operator and its perturbations occurring 
in a renormalization procedure, the following notion was introduced, see 
Definition 3.1 in [57] with L = P and P x = P 2 = L 

Definition 3.1. A symbol L(x,p; /i, h) will be called of strong type I if the 
following conditions are satisfied for all h G (0, h ) with some h > 0: 

(a) L depends analytically on p, G [—4,4] and takes real values for such p. 

(b) There exists e > such that 

(bl) L(x,p;p,h) is holomorphic in 

D £ = |(/i,x,p) GCxCxC:|/i|<4, \Qx\ + \%p\ < -|, 

(b2) for (p,x,p) G D £ , there holds 

\L(x,p; p, h) — (cosx + cosp + p)\ < e. 

(c) The following symmetry and periodicity conditions hold: 

(cl) L(x,p; p, h) = L(—p, x; p, h) = L(x, —p; p, h), 

(c2) L(x, p; p, h) = L(x + 27r, p; p, h) = L(x, p + 2tt; p, h) . 

By e(L) we will denote the infimum of e for which the above conditions hold. 
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In [251 ESI EZ] a detailed analysis was performed for pseudodifferential oper- 
ators associated with strong type I symbols. One of the results (appearing 
as an intermediate stronger result toward the proof of Theorem 0.1 in [27J) 

was 

Theorem 3.2. Let L(fi, h) be a strong type I symbol. There exists Eq > 
and C > such that if e(L) < Eq and if (27r) -1 /i is an irrational admitting a 
representation as a continuous fraction 

h _ 1 

2^ ~ 1 

ni H - 

n 2 H ■ 

n 3 + ... 

with integers nj satisfying nj > C , then the fi-spectrum of the associated 
family of operators Lh{fi) is a zero measure Cantor set. 

For operators H = (— zV + A) 2 + V with periodic potentials V, 

V(xi + 2ir, x 2 ) = V(xi, Xi + 2n) = V(xi,x 2 ), 

and constant magnetic fields, curl A = B, it was shown in several asymptotic 
regimes that the study of some parts of the spectrum reduces to a non- 
linear spectral problem of the above type for i? _1 -pseudodifferential operators 
with symbols close to V(x,p) for strong magnetic fields, see e.g. [28], or for 
-B-pseudodifferential operators with principal symbols coinciding with the 
first band function of the zero-field Hamiltonian (Peierls substitution) in 
the weak magnetic field limit, see [23 Appendix e], [23 p. 117], [2H], and 
earlier contributions by physicists mentioned e.g. in [8]. Hence, strong type 
I operators appear when considering potentials or first band functions close 
to /?(cosxi + cos x 2 ) + 7 for some reals (3 ^ 0, 7. 

Note that the expression ( 12.31) is a realization of the Krein resolvent formula 
for self-adjoint extensions, see e.g. [12], and can be viewed as a kind of 
the resolvent formula arising in the study of Grushin problems [33], which 
suggests a certain analogy with the case of regular periodic perturbations 
studied in [28]. Our aim here is to show that the spectral problem ( \2A\\ for 
periodic point perturbations can be also reduced in a certain regime (weak 
magnetic field and weak coupling) to the study of pseudodifferential operators 
with strong type I symbols. 



7 



4 Symbols associated with some discrete op- 
erators 



It is well known that the spectrum of the operator (13.21) as a set coincides 
with the spectrum of the discrete magnetic Laplacian acting on £ 2 (Z 2 ), see 
e.g. [2S1 pages 10-11], 

C h f{m, n) = e ihn f(m + 1, n) + e - lhn f(m - 1, n) + /(m, n - 1) + /(m, n + 1). 

The following theorem describes a similar correspondence for more general 
operators. This is essentially a suitable reformulation of the constructions of 
Section 6 in |28|. 



Theorem 4.1. Let Ch be a bounded linear operator on I 2 (7?) given by an 
infinite matrix (C(p, q)), p,q G Z 2 , with exponentially decreasing entries, 
\C(p,q)\ < ae~ b)[P ~ q]{ for some a,b > and all p,q G Z 2 , and satisfying 

C(p + k,q + k) =e- lhk2ipi - qi) C{p,q), p,q,keZ 2 , (4.1) 

with some h > 0. Then the spectrum of Ch coincides with the spectrum of 
the Weyl h- quantization of the (h-dependent) symbol T given by 



T{x,p)= c(m,n)e- imnh/2 e i{w, -" ! ' ] 



where c(m, n) = C ((0, 0), (m, n)) , m, n G Z. 

Proof. Let us start with preliminary constructions. Let T : M 2 — > C be a 
periodic smooth function, T(x,p + 2n) = T(x + 27r,p) = T(x,p). Consider 
its Fourier expansion, 



T{x,p) = t(rn,n)e t{mx+np \ 

i p2ty p2tt V / 

t( m ^ n ) = ——r / T(a;,p)e- i(ma;+np) dpcia;. 
(2vr)^ i Jo 

We will assume that the coefficients t(m, n) are exponentially decaying, i.e. 
that 

|i(m, n) I < ae~ bv/m2+ ™ 2 , for some a, 6 > . 
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Let us find an expression for the corresponding operator given by fl3.ll) . 
Let / £ <S(R). Substituting (g^D into (JHUD we obtain 



Thf(x) = -^- V t(m,n) / / e i(n P+ m( X+ y)/2] ^i P (x-y)/h f ( y)dpdy 

= — V t(m, n)e imi ' /2 / / e i{np+my/2) e ip{x - y)/h f(y)dpdy. (4.3) 
Note that for each i6l and n£Z there holds, in the sense of distributions, 

J e ip( x -y)/ h+ inp dp = _ 2n5 f^l + n ) . 

Jr h 
Hence Eq. ( 14.31) rewrites as 

T h f{x) = \ Yl t(m,n)e imx / 2 I e my ' 2 5{^ + n)f{y)dy 
= V t{m,n)e imx/2 (e imy/2 f{y)) 

* V / y=x+nh 

= t{m,n)e imx e imnh/2 f{x + nh), (4.4) 

and due to the exponential decay of the coefficients t(m, n) this extends by 
continuity to the whole of £ 2 (R). 

Let us return to the initial operator Ch- By the assumption (14. II) . 

C(p, q) = exp (ihp2(qi — Pi))c(q — p) for any p,q £ Z 2 . 

Hence 

Cfc/(P) = E ^ P2(9l - pi) c(g - p)f(q) = £ e lhp ^c(q)f(p + q). 

Therefore, commutes with the shift f(pi,P2) | — ► f{pi + 1,^2), and the 
Floquet-Bloch theory is applicable. 

Let us introduce the functions 

R 3 v ^ b n ((p) = ^2c(k,n)e ik(p , neZ. 
fcez 

All these functions are 27r-periodic and analytic in a complex neighborhood 
of R. Consider a family (indexed by 6 £ R) of operators C/j(0) acting in 
£ 2 (Z), 

Ch{0)g{m) = X^ b n (mh + 0)g(m + n), m £ Z . 

neZ 
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The operator- valued function 9 \— > Ch{9) is obviously continuous in the norm 
topology. Due to periodicity of the functions b n one has 

C h (9) = C h (9 + 2n) . 

Therefore, by the Floquet-Bloch theory, one has 

specC h = [J spec C h (9) . (4.5) 

0e[O,27r) 

Furthermore, for any 9 the operators Ch{9) and Ch{9 + h) are unitarily equiv- 
alent, 

C h (9 + h)=SC h (9)S- 1 , 
where S is the shift in £ 2 (Z), 

Sf(n) = f(n+1). 

This implies, cf. ( 14. 5p . 

specC h = [J spec C h (9) , 

6>S[0,h) 

which coincides with the spectrum of the following operator acting in 
£ 2 (Zx[0,ft)): 

Thu(m,9) = Ch{9)ug(m), UQ(m) = u(m,9), m G Z. 

Consider the map U : C 2 (R) -* £ 2 (Z x [0, h)) given by 

Uf(m, 9) = f(mh + 9) , (m, 9) e Z x [0, /i) . 

One can easily see that U is unitary. The operator := U~ l TjJJ is then an 
operator in £ 2 (M), 

f h f{x)=Y,K{x)f{x + nh)= c(m,n)e mx f(x + nh). 

nSZ m,nSZ 

Comparing this expression with (14.41) we arrive at the conclusion. □ 

Using Theorem 14.11 and the relations (12 .4p and (12.51) one arrives at the fol- 
lowing corollary. 
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Corollary 4.2. The set specify \ spec H h coincides with the z-spectrum of 
the family of h-pseudodifferential operators corresponding to the symbols 

L a {x,p;z,h) = Q{{0,0),{m,n);z)e- imnh / 2 e l(mx+np) + a' 1 (4.6) 



q(z, h) + a 1 + X(m, n; z, h)e 



i(mx+np) 



|m| + |n|>l 



where 



\(m,n;z, h) = F h ((0,0),(m,n);z) 

h{m 2 + n 2 ) \ /l z\ fl z h(m 2 + n 2 ) \ 

4 J \2 2h) \2 2h ' ' 2 /' K ' 
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and 

F h ( X ,y; z) = exp -*>(*« + «»> ) G „ (x , y , z) . (4 .8) 

The operators Lh, a associated with the above symbols L a can be viewed as 
effective Hamiltonians of the problem; we emphasize that these are given in 
an explicit form (the Fourier expansion). We are going to show that after 
suitable transformations one arrives at the study of strong type I symbols. 



5 Estimates for the Green function 

In the present section we will establish some estimates involving the Green 
function, the spectral parameter, and the strength of the magnetic field in 
the form which will be of importance below. The representation (12.1 j) works 
effectively near the Landau levels, see for example [20], while we intend to 
consider the negative halfline. So we prefer to use the corresponding heat 
kernel, i.e. the integral kernel of the semigroup e~ tHh , t > 0, see e.g. Section 
6.2.1.5 in [23], 



Ph(x v . t ) - h - exp ( - ih ^-y^ + y*) \ 



( h(x — y) 2 coth ht\ 
xexp(-^ ^ J. (5.1) 

To unify the notation, we denote by H the free Laplacian (i.e., the operator 
Hh with h — 0). The above expression ( 15.1 ft will be used for h > 0, and by 
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Pq we denote the heat kernel associated with the free Laplacian, 

1 / (x-y) 2 ^ 



P (x,y ] t) = —e* P (- u 
For the Green functions we will use the integral representation 

POO 

G h (x,y;z) = e zt P h (x,y;t)dt, h > 0, (5.2) 
Jo 

which is valid at least for $lz < h. For the function Fh introduced in ( 14.8ft 
one has 



oo 

zt 



F h (x,y;z)= e zt \Ph(x } y-t)\dt, h > 0. (5.3) 
J o 

Recall that the Green function Go of the free Laplacian is known explicitly 
as well, 

G (x,y;z) = —K (y^z\x-y\), (5.4) 



2tt 

where K Q is the modified Bessel function of order zero. Here and below we 
fix the branch of the square root which is a continuation of the usual square 
root on the positive half-line. 

We will also use the integral kernel G^\x,y; z) of the operator (Hh — z)~ 2 , 
h > 0, ^tz < 0, and the function 

ff >(*, y; z) = exp (^l^lllA) G f (x , y , z) . 

By the Hilbert resolvent identity, there holds 

G h 2 \x,y;z), (5.5) 



dG h (x,y;z) r<t2 . 



dz 

and, at the same time, 

a -^^ = FP( X ,y;z). (5.6) 

Using the identity K' (w) = -K x {w), see Eq. (9.6.27) in PQ, and Eq. (l5~4"j) 
above, we obtain 

G { *\x, y- z) = K X (V=i \x-y\). (5.7) 

On the other hand, using 

POO 

(H h - z)- 2 = / te- t{Hh - z) dt, 3^ < 0, h>0 
Jo 
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we arrive at 

/■oo 

Gf\x,y;z)= / te zt P h (x,y;t) dt, ^z<h, h > 0, (5.8) 

/■oo 

F ( ^\x,y]z)= I te zt \P h (x,y;t)\dt, Viz < h, h > 0. (5.9) 

The following two-side estimate will be of importance. 
Lemma 5.1. For h > there holds 

exp ( - ht - k{ - X ~ y) )P (x,y;t) < |P fc (x,2/;t)| <P (x,y;t) 

/or any t > and i,!/6l 2 . 

Proof. The right-hand side inequality is nothing but the Kato (diamagnetic) 
inequality valid for much more general magnetic fields, see e.g. §2 in [UJ, so we 
only need to prove the left-hand side inequality. We will use the elementary 
estimates t 

> 1 and — t —— > 1 for t > 0. (5.10) 

Using (15.1 Oft one has 

h e~ ht 2ht e~ ht 

> 



sinh ht t 1 - e" 2/l * ~ t 
At the same time, using ( 15.101) again, one obtains 

, /, 2 \ , 1 2ht ,1 

hcothht = h[l + J = h H ^- < h+ -. 

V e 2ht - 1/ t e 2ht - 1 ~ t 

Therefore, 

, . , , h ( h(x — y) 2 coth 

>^expf-/.t- /l(X -^ )2 -^° 
= exp ( _ ^ _ j p o(x? y . f y 

Using (15 .2p and (15.81) one arrives at the following corollary. 



□ 



13 



Corollary 5.2. For x, y G M 2 , x 7^ y, and real z, z < h, one has 



\G h (x, y;z) \ > G (x, y;z-h) exp 



h(x - yf 



and, for any x,y G M 2 and real z with z < h, there holds 



\G h (x,y;z)\ > (x,y; z - h) exp - 
Lemma 5.3. For h > and ?R,z < t/iere holds 
\F h (x,y;z) -G (x,y;z)\ 



h(x — yf 
4 



< Go (a;, ?/; Jfc) — Go(:r, 1/; 9fo — h) exp f — 



h(x - y f 



(5.11) 



and 



Ft 2 \x,y;z)-Gf ) (x,y;z)\ 



(2), 



< Gj, 2 ' (x, y; 3fcz) - G z; (x, y ; 3fc - h) exp ( - 
Proof. Using ( 15. 2p and ( 15.31) one has 



/i(x — y) s 



(5.12) 



l-Ffcfol/jz) - G (x,y; z)\ 



POO 

/ e z *(|P ft (x !2 /;t)| - P {x,y;t))dt 
Jo 

Ph(x,y,t) \ - P (x,y;t) 



< I \< 
'0 



/i(x — yf 



P (x,y;t). 



Estimate | e zt | = e 3 * 2 * and, using Lemma 15.11 

|-Ph(aM/;t)| - Po(x,y;t) = P (x,y;t) - \P h (x,y;t) 

< P (x,y;t) - exp ( - /it 
We obtain, using ( 15. 2ft again, 
|F h (a:,j/;z) - G? (a;,y;«)| 

< J°° e ** (P (x, y; t) - exp ( - ht - k{x ' ~ y)2 ) P (x, y; t))dt 

= J™ e^P (x, y; t)dt - exp ( - J™ e (^)tp o{Xj y , t ) dt 



G (x, y; $tz) — Gq(x, y; $lz — h) exp ( — 



h(x - yf 



Eq. ( 15.12ft is obtained in the same way using ( 15.81) and ( 15. 9ft . 



□ 
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In what follows we need the function 



q (z) := lim (g o (x, y;z)- — log ■= — - 
x-+y \ 2ir \x — 



y\ 



(5.13) 



Using Eq. (9.6.13) in [T], 



lim (K (wr) + logr) = — logw; + log2 + ip(l), u>eC\(— oo,0], 

r— >0+ 



we obtain 



log(-z) -log4-2V>(l) 

47T 



(5.14) 



Everywhere we take the principal branch of the logarithm, i.e. the holomor- 
phic extension to C \ (— oo, 0] of the usual logarithm on (0, +oo). 

Lemma 5.4. For h > and tflz < there holds 



\q(z,h) - q (z)\ < ^- log ■ ~ 



dq(z,h) dq (z) 



dz 



dz 



< 



$lz 
h 



4tt3?z(3?z - h) 



(5.15) 
(5.16) 



Proof. Consider first (15.151) . Due to the logarithmic on-diagonal singularity 
of the Green functions one has 



lim ( G h (x, y; z) - F h (x, y; z) 

x—>y \ 



Ito (! _ exp ( j*^zJ^±g) )) Gk (z,y; z) = 0. 



Hence, using (12.21) . 

q(z, h) = lim (F h (x, y;z)- — log , — - — -) . 

x-*y \ Its: \x — y\) 

Substituting now (15.171) and (15.131) into (15. lip we see 

\q{z,h) -qo(z) \ = lim F h (x,y;z) -G (x,y;z) 



(5.17) 



x^y 



( h(x — y\ \ 
< lim ( Go (x, y; $lz) — Go (x, y;^Rz — h) exp ( ) ) 

x^y V 4 / 



lim (g (x, y; %tz) - ^- log - — - 
x^y \ 2ir \x — 



y\ 



- lim (g (x, y;^z - h) - — log j- ^ 

1 / / h(x - y) 2 N \ , 1 
+ hm — 1 - exp ( ) log : 

x-*y 2n \ v 4 ' ) \x — y\ 



f h(x-y) 2 , 
exp ( ) 

g (Sfcs) - qoi^z - h), 
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and now one can use the explicit expression (15.141) . 

As for (15.161) . the expressions (12.21) for q(z,h) and (15.131) for qo(z) together 
with the analyticity of Gh in z (h > 0) imply 

= &h \ X i X 'i Z ) = \ X i X \ Z )i = \x, X] z), 

where x G M? is arbitrary. Therefore, using ( 15.121) . one obtains 



dq(z,h) dq (z) 



dz dz 
and, by ( 15.141) . one has 



< G^ (x, x; fUz) - Gf (x, x; $lz - h) 
dq (X) 1 



dqo(X) 



dX 



X=?Rz 



A=J?z-/i 



4ttA 



which finishes the proof. 



□ 



Below we will use intensively the following well-known asymptotics, Eq. 
(9.7.2) in P: 



K (w) = J^e- w (l + O(w- 1 ))., 



7T 



— > oo, | argu?| < — . (5.1£ 



Let us prove first several estimates from above for the Green function and 
its derivative. 

Lemma 5.5. For any d, S > there exists E > such that for \x — y\ > d, 
3lz < —E, and h > there holds 



\G h (x,y;z)\ < 



7T 



1 + 8 J 

2tt \l-4Mz\x-y\ 2 



-Uz\x-y\ 



Proof. Using ( 15.21) and the right-hand side inequality of Lemma 15. 1[ for 
3lz < we obtain 



\G h {x,y,z) 



POO 

/ e zt P h (x,y;t)dt 
Jo 



poo poo 

< / \e zt \ ■ \P h {x,y;t)\dt< / e* zt P {x, y; t)dt = G {x, y; toz). 
Jo Jo 

Substituting here ( 15.41) one obtains 



\G h (x,y;z)\ < — K {\T^z \x - y\), 

1 1 Z7T 



and the asymptotics ( 15. 18ft gives the result. 



□ 
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Lemma 5.6. For any d, 5 > there exists E > such that for \x — y\ > d, 
IRz < —E, and h > there holds 



dG h {x,y;z) (1 + 5) Jtc 2 \x - y\' 



dz 



< 



Air 



-4(3?z) 



exp ( — \J — $tz \x — y\). 



Proof. Using Eqs. (E3J), (15.81) . and the right-hand side inequality of 
Lemma 15.11 for 9ftz < we obtain 



dG h (x,y;z) 



dz 



G^{x,y;z)\ < I t\e zt \ ■ \P h (x,y;t)\dt 



POO 

< / te !ltzt P {x,y;t)dt = G i o ) {x,y;$lz). 
Jo 



Hence, using (15.71) we arrive at 

dG h (x,y;z) \x-y\ 



dz 



< 



AttV^Wz 



K^y/^lx-yl), 



and it is sufficient to use the asymptotics 

Ki(w) = J — e~ w (l + Oiw' 1 )), \w\-*oo, |argw|< 
V 2w 

see Eq. (9.7.2) in p. 



□ 



Let us pass to estimates from below, which are of crucial importance. 

Lemma 5.7. For any x, y E M 2 , x ^ y, and any c G (0,1) there exist 
E, 7, ho > (which depend on x, y, c) such that 



\G h (x,y;z)\ > 



1 - c 



7T~ 



-®z\x-y\ 



2tt y -4$tz\x-y\ 2 
for $lz < —E, \^sz\ < 7, and h 6 (0, ho). 



Proof. By Lemma 15. 3[ 

\G h (x,y;z) \ = \F h (x,y; z)\ 

> \G (x,y;z)\ - \F h (x,y;z) -G {x,y;z)\ 

> \Gq(x, y;z) \ + G (x, y; — h) exp 



\Go(x,y;z) \ + G (x,y;$lz - h) (^exp 
+ (G Q (x, y; ^kz — h) — G (x, y; $lz)). 



h(x - y) : 
4 

h(x - y f 



- Go{x,y\ 3fcz) 
(5.19) 
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Due to (I5.18P and Lemma l5~5l for any 5 £ (0, 1) one can choose E(S, x, y) > 
such that 



1-5 



2tt 



7T 



4z|a; — y\ 2 
1 + 5 



< 



7T" 



2tt y -4Kz|x-y| 2 

The constant 5 will be chosen later. 
One can find ho = h (5, x, y) such that 



< \G (x,y; z)\ 

e -^W-z\x- y \^ sa z< -E(6,x,y). (5.20) 



ho(x - yf 
vxp \ ) - 1 



then, by the right-hand side inequality in ( 15.20D . 



G (x, y; Kz-h)l exp ( - ^ , ^ ) I 



<5, 



(5.21) 



4 



< 



7T 



5(1+5) 

2tt y 4(/i-^)|x-y| 2 
5(1+5) 



< 



7T 



-Vh-9iz\x-y\ 



-$tz\x-y\ 



2tt \/-4^|x-y| 2 
for 3fb < -E(5,x,y), h e (0,h (S,x,y)). (5.22) 



At the same time one has 

exp ( — \f^z\x — y\) exp ( — tft^— z\x — y\) 



> 



exp ( 



-z\ \x 



- y\) _ exp ( 



z\ \x 



-y\) 



z\ 

> 



hence for \^sz\ < 7 one can estimate 



^sz\ 



exp 



(-y/\®z\ + \Qz\\x-y\), 



ex.p(-y^z\x-y\) > ===== exp (- y^T\^z]\x-y\). (5.23) 

V7+ 1^1 
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For t > E > and 7 > one can write 



VtT^-Vi=—=L t<^t< " 



Hence, for any a, 7 > and t > E > 0, one has 

exp ( — < exp ^ X_ J exp ( — oty/t + 7) . (5-24) 

Analogously one has for t > E > and 7 > 0, 



^7 + *=^/l + j^< (l + j) v^< + (5-25) 



Hence, assuming that 7 = j(5, x,y) satisfies 

l\x-y\ .... 7 . r 

exp < 1 + <), — < d, 

2v^ ~ E ~ 



and substituting (15.231) . (I5.24p . and (15.251) into the left-hand side inequality 
in (I5.20p one arrives at 



7T 2 



\GJx vz)\> — — 4 4 / - e -V=*5\x-y\ 

for 3lz < —E(5,x,y) and < j(8,x,y). (5.26) 
The estimates (15.201) give also, for Jfc < — E, 
\Gq(x, y; $tz — h) — Gq(x, y; diz) 



< 1 + 5 </ ^ c ->/=55|s-»| 



2tt V -43fc |ac - y\ 2 



1 — 5 J 71 



2 



2tt V A(h - 9fcz) \x - y\ 2 



e -Vh^\x-y\ 



" 2tt\V -4$tz\x-y\ 2 V 4(/i - Uz) \x - y\ 2 



+ 



5 J 71 



2 



4 11 -V^\x-y\ 

7i \I -4Mz\x-y\ 2 
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Expand 



^2 2 



4/ ->/=IE[x-v| _ 4/ « ^Sil^-yl 

-43?z|x-y| 2 V 4(/i - 3fcz) |x - y| 2 



7T 2 



-43?z|x-y| 2 



e -V=55|a:-i/| _ e -Vh=^\x-y\ 



, | ,■*/ 71-2 .*/ ![! \ -VS=S51*-y| 

1 l; -4ftz|a;-y| 2 y 4(fc - 3&) |x - j/| 2 / 

Using (15.241) we estimate, for h G (0, ho) and 5fc < — E, 



exp ( — V— 3lz\x — ?/|) — exp ( — yh — ?R,z\x — y\) 

1 J exp ( — \Jh — ^lz\x — y\) 



^ i h\ x - y\ 

< I exp 



2VE 

h \x - y\ 



< ^exp ° X 1 j exp ( — y/— 3t.z\x — y\). 



2VE 

At the same time one has, for h G (0, ho) and 3lz < —E, 



7T 2 



-A$lz\x-y\ 2 V 4(/i - $lz) \x - y\" 



u h-?ftz \ 4 / 7T 2 



-3f?2 J V 4(/i-8fcz) |z-y| 



2 



< ____ 4/ /iQ 4/ TT" 



|3fcz| y 4|3te| \x - y\ 2 ~ E y 4|9fcz| |x - y| 2 ' 
Let us choose /io = h (5,x,y) in such a way that (15.2ip still holds and that 

fh \x-y\\ h . 

exp — — — 1 < 0, — < 0, 

V 2v ^ / E 

then 



I G (x, y; ttz-h)- G (x, y; Sfc) | < ^ _ y|2 e^l^l 

for 3fc < -E(5,x,y), \$Sz\ < j(6,x,y), h G (0, h (S,x,y)). (5.27) 
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Substituting now the estimates (I5.22p . (I5.26p . and (I5.27P into (15.191) . one 
shows that 



e -V=S5|*-»| 



-4!Rz \x — y\ 2 

for < —E(S,x,y), \$iz\ < j(5,x,y), and h e (0, /io(5, 2/)), which yields 
the result due to the arbitrariness in the choice of 6. □ 



6 Estimates for the effective Hamiltonians 

In this section, we are going to analyze the symbols L a from Corollary 14.21 
One can rewrite (14.61) as follows: 



L a (x,p; z, h) = q(z, h) + a 1 + a(z, h)( cosx + cosp) + W(x,p; z, h), 

a{z,h) := 2X(1,0; z,h), 

W(x, p;z,h):= A ( m ' n 5 z ' h ) e l{mx+np) . 

m,n£Z, 
|m| + |n|>2 

First note that due to the integral representation, Eq. (13.2.5) in p], 

POO 

T(a)U{a,qx) = / e-* t t°- 1 (l + t) c -°- 1 cft, > 0, 
Jo 

and to the formula (12. ip for the Green function one has 

1 / ih(x 1 -y 1 )(x 2 +y 2 ) h(x-y) 2 * 



G h (x,y; z) = — exp ( 



x / e 2 



Therefore, by (14.71) . there exists C = C{h) > such that for Jfc < one has 

I A(m, n; z, h)\ < C exp ^ J for all m and n, \m\ + |n| > 0. 

Hence L a is an analytic function of (x,p) G C 2 . 
Note that due to Corollary 15.21 and (15.41) one has 

i-h/4 



\a(z,h)\>— — K (Vh — z), h > 0, zE(—oo,h). 
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As the function K (w) is strictly positive for positive w, see e.g. Section 9.6.1 
in [T], one can conclude that a(z,h) is non-zero for all positive h and real 
negative z, and Corollary 14.21 gives 

Lemma 6.1. The negative spectrum of Hh a coincides with the negative z- 
spectrum of a family of h-pseudodifferential operators associated with the 
symbols 

M a (x,p; z, h) = 1 L a (x,p; z, h) 
a{z, h) 



with 



= m a (z, h) + cos a; + cosp + N(x,p; z, h) (6.1) 

q(z,h) + a~ l 
= a (z,h) ' 

N{x,p;z,h)= T A(m ; n ;^ } e^rt. ^ 

|m| + |n|>2 

Lemma 6.2. Ta&e /io > and R > 0. TTien one can /md Eq > wift 
t/ie following property: for any E > E there exists otE > stzc/i t/iat /or 
a G (0, cue), h £ (0, ho), 3fc < —E the condition | g(z, /i) H-a" 1 1 < R\a(z, h) | 
guarantees that $tz < —E and \^sz\ < Ah . 

Proof. Take an arbitrary E' > 0. Due to Lemma 15. 5\ there exists C such 
that, for §lz < —E' and h > 0, we have the estimate 



\a(z,h)\ < C'exp ( - V^Rz) . 
Therefore, the condition 

\q(z, h) + a~ 1 \ < R \a(z, h)\ 

implies 



\q(z,h) + a~ 1 \ < Cexp{-V-$tz), C = C'R. 
Using the triangular inequality, 

\qo(z) + a" 1 ) < \q(z, h) + aT 1 ) + \q(z, h) - q (z)\, 

and the estimate ( 15. 15ft we obtain 

|»M + »-| <l log ^ + Ce~^ < ^ + (VVH. 
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Substituting here the explicit expression (I5.14p we obtain 



log \z\ + i arg(— z) — log 4 — 2-0(1) 

An 



+ a' 1 



< 



h 



An\^z\ 



In particular, considering the real and the imaginary parts separately, 



log \z\ - log 4 - 2^(1) 



An 

axg(-z)| 



< 



+ a' 
h 



< 



h 



An ~ An\^z\ 
Using (16 .4p . let us choose Eq > E' such that 



An\Rz\ 



Ce 



-»2 



(6.3) 
(6.4) 



n 



| arg(-z) | < -, $iz < -E , h £ (0, h ) 
log2£ - log 4 - 2i/>(l) > 0. 
Then, under the same conditions on z and h, one has 

|9fcz| < |z| < 2|3fcz| , 

and from (16.31) one can get 



h 



An\^z\ 



+ Ce 



^ | log2|^|-log4-2^(l) 
An 



(6.5) 



h 



With 
B:=swp\- 

the estimate (I6.5P gives 



2-0(1) + log4 

47T 



: 3fc < -£ , /i £ (0,/i )} 



c^ 1 < 5 + 



log2|9Kz| 

47T 



(6.6) 



Therefore, for any E > Eq and a £ (0, a^) with 

47T 



47rS + log 2E 
the condition ( 16.61) gives $lz < —E. 

Using ( 16. 4p again, one can (if necessary) increase £7o to obtain 

2h 



arg(— z)\ < 



\Mz\ 
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^Z < -En 



and 

[ 2/i 

tant < 2t for t E 0, — - 
V E 

then one has 

\Qz\ = | tanarg(— z) \ ■ \$tz\ < 4/i 
for <-Eq. □ 

Lemma 6.3. Take ho > and R > and Zet £"o &e t/ie corresponding 
constant from Lemma lUTB . There exists a>o > suc/i £/&a£ 7 for any h G (0, ho) 
and a e (0, cto) > equation 

q(z, h) + a -1 = /ia(z, h) (6.7) 

wjf/i respect to z, §lz < —Eq, has a unique solution z = ( a (fi, h) holomorphic 
in \x for complex \x with \fi\ < R, and this solution is real for real \i. 

Proof. For complex ft and a > denote by z(ft,a) the solution z to the 
equation 

q (z) + a" 1 = ft. 

Due to the explicit formula (15.131) for q one has an explicit expression 

z{ft, a) = - exp (bg 4 + 2^(1) + - 4tt/3J ; (6.8) 

in particular, 

3te(/3,ot) = - exp ^log4 + 2^(l) + ^ - 4tt3?/5J cos(4ttS/3). 

Now note that if /3 = /3(a, /i, h) satisfies 

ft = q (z(ft,a)) - q(z(ft,a),h) + fxa(z(ft,a),h), 

then C«(/ i , h) := ,z(/3(a, /i, h), a) is a requested solution for (16. 7\j . 

To show the existence of a unique ft (a, n, h) we use the fixed point theorem. 
We are going to analyze the map 

ft h-> /i) := qo(z(ft, a)) - q(z(ft, a),h)) + /ia[z(ft, a),ti), 

and to show that it is a contraction for \ft\ < 1/16 if a is sufficiently small. 
Namely, we will prove that one can choose ao in such a way that 

|tf(/3,a,/i,/i)| < ±- for \ft\ < — , ae (0,a ), < R, h E (0, /io), (6.9) 
1 1 32 16 
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and that 



a, A*, /i) - ^(/3 2 , a, A 4 , h) 



< 



|/5i-/? 2 | 



for |/3 2 | < — , a G (0, a ), < R, h £ (0, /i ), (6.10) 
16 

then, by the fixed point theorem, the equation (3 = \I/(/3, a, ji, h) has a unique 
solution with \/3\ < 1/16 depending holomorphically on //, which can be 
obtained throught by iteration /3 n+ \ = ^(/3 n ,a,fi,h) with, say, (3q = 0, and 
taking the limit j3 = lim/5 n . If /i is real, then all f3 n are real as well, hence so 
is (3. 

Note that for \(3\ < 1/16 one has 



1 4-7T 

9fcz(/3,a) < --^ exp (log 4 + 2^(1) + — - 4vr 



(6.11) 



To show (16.91) we note that, by Lemma I5T41 and Lemma 1531 there exists E\ > 
such that for $tz({3, a) < —E\ and < R one has 



|#(/3,a,/i,/i)| < 



h 



+ 



R 



exp 



(-y/\Xz(P,a)\). 



4vr|^(/3,a)| |^(/5,a)|V4 
Therefore, there exists E 2 > E\ such that 

\^((3,a,fM,h)\ < — for ft, G (0,/i ), H < i?, 3te(/3, a) < -E 2 , 

and it remains to use (16. lip to choose a% > in such a way that 9fc(/3, a) < 
-E 2 for a G (0,ai) and < 1/16. 

Let us prove now the estimate ( 16. 10(1 . One has 

a, - *(/3 2 , a, //, ft) < $(a,/i,/i)|/?i -/? 2 |, 

<9^(/3, h) 



where 



$(a, /i, h) := sup 

|/3|<1/16 

Hence it is sufficient to show that 

d^((3,a,n, h) 



d(3 



d/3 



< 1/2 



(6.12) 



for \f3\ < 1/16 and small a. 



25 



We note first that, by 
such that 



and (RETT]) , for \p\ < 1/16 there exist B,C > 



dz(P, a) 



dp 



An i 
< Bexp— < C\Uz(p,a) 
a 1 



(6.13) 



By Lemma EH and Lemma l5\6l there exists E3 > such that for 3lz(P, a) < 
—E 3 one has 



dp 



< 



dz(f3, a) 



dp 



h 



4tt mz(P,a) 



+ 



R 



1 3/4 



< 



C7i 



47r 3fc(/?, a 



T + C J R|^(/3,a)| 1/ Vv^(^, 



where the last estimate is due to ( 16.131) . Hence one can find E4 > E 3 such 
that (1632]) holds provided Uz(P, a) < -E 4 for \p\ < 1/16 and h G (0,h ), 
and one can again pick a 2 G (0, «i) in such a way that 5fc(/3, a) < — E4 for 
< 1/16 for a G (0,0^2) using the estimate (16. lip . 

Now it remains to show that for sufficiently small a any solution z to (16 .7p 
really satisfies \qo{z) + a -1 < 1/16; this will mean that the above solution 



( a is unique. For any such solution z one has \q(z, h) 



a 



-11 



< R\a(z,h)\. 



Combining this with the triangular inequality and the estimate of Lemma 
15.41 one arrives at 



\qo(z) + 01 1 \ < \q(z, h) + a x \ + \q(z, h) - qo(z)\ < R \a(z, h)\ + 



h 



4tt|3?z| 



By Lemma I6.2[ for any E > E we can find G (0, 02) such that for a G 
(0,a£) and h G (0, /i ) the solution satisfies Sftz < — E. Take sufficiently 
large to have 



R\a(z, h)\ 



he 



4n\$tz\ 



< 1/16 for 3fc < -E, 



which is possible due to the estimate of Lemma l5\6l for a(z, h). Denoting the 
corresponding value of by a , this completes the proof. □ 

Below we will need the fact that a(z, h) does not vanish in a certain tubular 
complex domain, and not only on the negative halfline. Namely, it follows 
from Lemma [5. 71 that one can find h* > and E* > such that : 



a(z, h) ^ , if 3ftz < -E* , \$Sz\ < Ah* and h G (0, K 
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Note that in this domain the equation (16.71) reads simply as m a (z, h) = /i. 

The following theorem summarizes all previous considerations and shows the 
applicability of the Helffer-Sjostrand methods [231 l2"B"t 127]. 

Theorem 6.4. For any E > there exist «o > and ho > such that for 
a G (0, a ) and h G (0, h ) one has 

specify n (-oo,-E ) = {z < -E : m a (z,h) G //-spec P hj(X (//)}, 

where m a is given in (|6.2p . P^ Q zs a pseudodifferential operator with a strong 
type I symbol P a (x,p; /i, h) given by Eq. (16.141) below, and lim^Q e(P a ) = 0. 

Proof. By Corollary 14.21 and Eq. ( 16. ip . the negative spectrum of Hh, a co- 
incides with the z-spectrum for the symbols M a (x,p; z,h). Note that the 
symmetry conditions (c) in Definition 13.11 are satisfied due to the specific 
form of the Fourier coefficients, see (14. 7p . so we will be concerned below with 
the conditions (a) and (b) of Definition 13.11 

Due to Lemma [531 and the estimate (15.151) . for any segment [a, b] C (— oo, 0) 
there exists a constant C > such that h)\\ <C for z G [a, b] and any 

h. By ( 12. 4p . this means that Hh, a has no spectrum in [a, b] for |a| < 1/C. 
Therefore, we can assume without loss of generality that Eq satisfies the 
assumptions of Lemma 16.21 for R = 5 and h = h* , and that E > E*. 

Using Lemma 16.31 one can choose a\ > such that for h G (0, h*) and 
a G (0, «i) the equation m a (z, h) = \i with respect to z with 9ftz < — E has 
a unique solution z = (ait 1 , h) for all \x with < 4, and this solution is a 
holomorphic function of /i and is real for real [i. 

Let us consider the symbols 



P a (x,p;p:,h) := M a (x,p;( a (fi,h),h). 



(6.14) 



Clearly, 



P a (x,p; fi, h) := /i + cos a; + cosp + T a (x,p; /i, h), 



where 



T a (x,p;n,h):= ]T t a (m,n; ^, h)e l ^ mx+n ^ 



m,n£Z, 
|m| + |n|>2 



and 



t a (m, n; fi, h) : 



\(m,n;( a ((jL, h),h) 
a(( a (fi, h),h) 
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As previously noted, the symbols M a and hence T a are analytic with respect 
to x and p in the whole complex space. Furthermore, for real p the solu- 
tions ( a (p>, h) are also real, hence, by (14.7p . the coeffcieints t a (m,n; p, h) are 
real as well and, as noted above, satisfy t a (m,n; p,h) = t a (—m,—n;p,h). 
Therefore, P a takes only real values for real p. 

Now take an arbitrary e > 0. By Lemmas 15.51 and 15.71 there exists E\ > Eq 
and h G (0, h*) such that 

\t a (m, n; p, h) | < 2 exp ^ — a/— tftCaip, h) (V m 2 + n 2 — l) j . 

for 3t( a (p, h) < —Ei, 9?£ a (jtf, ^) < 4/io, ^ G (0, /^o), and \m\ + |n| > 1. Then 
one can estimate for \^sx\ + \^sp\ < s^ 1 (hence for \^sx\ < s^ 1 and \Qp\ < s^ 1 ): 

\T a (x,p; p, h)\ < 2 ^2 exp + — - y/-^( a (p, h) (y/m 2 + n 2 - l) 

m,n£Z. 
\m\ + \n\>2 

Therefore, there exists E > E x such that 

\T a (x,p; p,h) \ < e for \$$x\ + < - 

provided 

< 4, $l( a (p,h) < -E, \Q( a (p, h)\ < 4h , he(0,h ). 

Now using Lemma [6T21 we can choose «o € (0, a\) such that ?R( a (p, h) < —E 
and |^ Q (/x, /i) | < Ah for all a G (0, cko), ^ £ (0,/io), and < 4. This 
means that for all a G (0, ato) and /i G (0, /io) the symbol P a (p, h) is of strong 
type I with e(P a ) < s. □ 

Using Theorem 13.21 one can get the following corollary. 

Corollary 6.5. For any E > there exist a > and C > suc/i t/iat /or 
a G (0, «o) and sufficiently small irrational (27r) _1 /i admitting a continuous 
fraction expansion 

h_ 1 

2^ ~ 1 

ni H j 

n 2 H ■ 

n 3 + ... 

with integers rij satisfying rij > C, £/ie spectrum of Hh, a in ( — oo, — -Eo) ^ a 
zero measure Cantor set. 
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Proof. For sufficiently small e(P a ) and h as described, the set //-spec Ph^in) 
is a zero measure Cantor set in virtue of Theorem 13.21 In view of Theorem 
16.41 it remains to emphasize that z \— > m a (z, h) is an analytic topological iso- 
morphism for sufficiently small a and h (see Lemma 16. 3p , hence the preimage 
of a zero-measure Cantor set is a zero-measure Cantor as well. □ 

Remark 6.6. In [TT] another solvable model was considered, namely, peri- 
odic quantum graphs with magnetic fields, and the expression for the spec- 
trum obtained there in Theorem 13 is very close to the one of Theorem 16.41 
Nevertheless, the essential difference was that the term T a in the construc- 
tions analogous to Theorem 16.41 would just vanish, hence one had a reduc- 
tion to the usual spectral problem for the discrete magnetic Laplacian. It 
is worthwhile to note that, even with such a simplification, the spectrum of 
non-isotropic quantum graphs was studied only numerically [21 J up to now. 
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